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Abstract 

In this article, we investigate the contribution of the higher-twist Feynman diagrams to the large-py 
inclusive single meson production cross section in photon-photon collisions and present the general 
formulas for the higher-twist differential cross sections in case of the running coupling and frozen 
coupling approaches. The structure of infrared renormalon singularities of the higher-twist subprocess 
cross section and the resummed expression (the Borel sum) for it are found. We compared the 
resummed higher-twist cross sections with the ones obtained in the framework of the frozen coupling 
approach and leading-twist cross section. We obtain, that ratio R = for all values 

of the transverse momentum px of the meson identically equivalent to ratio r={A^^Y'^^ / (A^.^)^ . 
It is shown that the resummed result depends on the choice of the meson wave functions used in 
calculation. Phenomenological effects of the obtained results are discussed. 
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I. INTRODUCTION 



Exclusive processes involving large momentum transfer are among the most interesting and 
challenging test of quantum chromodynamics (QCD). The framework for analyzing such pro- 
cesses within the context of perturbative QCD (pQCD) has been developed by Brodsky and Lep- 
age [1,2], Efremov and Radyshkin [3], and Duncan and Mueller [4]. They have demonstrated, to 
all orders in perturbation theory, that exclusive amplitudes involving large momentum transfer 
factorize into a convolution of a process-independent and perturbatively incalculable distribu- 
tion amplitude, one for each hadron involved in the amplitude, with a process-dependent and 
perturbatively calculable hard-scattering amplitude. 

The hadronic wave function in terms of quark and gluon degrees of freedoms plays an 
important role in QCD process predictions. For example, knowledge of the wave function allows 
to calculate distribution amplitudes and structure functions or conversely these processes can 
give phenomenological restrictions on the wave functions. 

During the last few years, a great deal of progress has been made in the investigation of the 
properties of hadronic wave functions [5- 17]. 

The standard approach to distribution amplitudes, which is due to Brodsky and Lepage[14], 
considers the hadron's parton decomposition in the infinite momentum frame. A conceptually 
different, but mathematically equivalent formalism is the light-cone quantization[15]. Either 
way, power-suppressed contributions to exclusive processes in QCD, which are commonly re- 
ferred to as higher-twist corrections. The higher-twist approximation describes the multiple 
scattering of a parton as power corrections to the leading-twist cross section. 

Among the fundamental predictions of QCD are asymptotic scaling laws for large-angle 
exclusive processes [18-22]. QCD counting rules were formalized in Refs. [19,20]. 

The frozen coupling constant approach can be applied for investigation, not only exclusive 
processes, but also for the calculation of higher-twist contributions to some inclusive processes, 
for example as large -pr meson photoproduction [23], two-jet -l-meson production in the electron- 
positron annihilation [24]. In the works [24,25] for calculation of integrals, such as 

j~/"-w;)^<;-^''d. (1.1) 

the frozen couphng constant approach was used. According to Ref.[25] should be noted that 
in pQCD calculations the argument of the running coupling constant in both, the renormaliza- 
tion and factorization scale should be taken equal to the square of the momentum transfer 
of a hard gluon in a corresponding Feynman diagram. But defined in this way, as{Q'^) suffers 



from infrared singularities. For example in our work [26], equals to X2S and —X\u, where s, 
u are the subprocess's Mandelstam invariants. Therefore, in the soft regions a; — )■ 0, integrals 
(1.1) diverge and for their calculation some regularization methods of q;s((5^) in these regions 
are needed. In Ref.[27], the authors investigated the phenomenology of infrared renormalons in 
inclusive processes. The dispersive approach has been devised to extend properly modified per- 
turbation theory calculations towards the low-energy region [28]. Connections between power 
corrections for the three Deep Inelastic Scattering sum rules have also been explored in [29]. 

Investigation of the infrared renormalon effects in various inclusive and exclusive processes is 
one of the most important and interesting problems in the perturbative QCD. As we know the 
word "renormalon" first appeared in Ref.[30]. A singularity in the Borel parameter- is called 
a renormalon. It is known that infrared renormalons are responsible for factorial growth of 
coefficients in perturbative series for the physical quantities. But, these divergent series can 
be resummed by means of the Borel transformation [30] and the principal value prescription 
[31], and effects of infrared renormalons can be taken into account by a scale-setting procedure 
as(<5^) (^s{^xp{f {Q'^))Q'^) at the one-loop order results. Technically, all-order resummation 
of infrared renormalons corresponds to the calculation of the one-loop Feynman diagrams with 
the running coupling constant as{—k^) at the vertices or, alternatively, to calculation of the 
same diagrams with nonzero gluon mass. Studies of infrared renormalon problems have also 
opened new prospects for evaluation of power- suppressed corrections to processes characteristics 
[32]. Power corrections can also be obtained by means of the Landau-pole free expression 
for the QCD coupling constant. The most simple and elaborated variant of the dispersive 
approach, the Shirkov and Solovtsov analytic perturbation theory, was formulated in Ref.[33]. 
The kx factorization theorem has been widely applied to inclusive and exclusive processes in 
perturbative QCD. This theorem holds for simple processes, such as deeply inelastic scattering 
(DIS) and Drell-Yan production [34]. 

A full twist 3 treatment of p-electroproduction in kr -factorisation is possible[35]. It relies 
on the computation of the 7^ — pr impact factor at twist 3 including consistently all twist 
3 contributions, i.e. 2-parton and 3-parton correlators. This gives a gauge invariant impact 
factor, and an amplitude which is free of end-point singularities due to the presence of kx ■ 

An additional, general property is that the singularities in B[u\ occur at integer-and some- 
times half integer- values of u. This corresponds to the fact that alternative definitions of the 
sum of the series differ by integer-or half integer-power of A^/Q^. These ambiguities must be 
cancelled by nonperturbative power corrections, and they can therefore serve as a perturba- 
tive probe of such effects. In the absence of an operator product expansion, the renormalon 



technique often provides a unique window into the nonperturbative regime: by identifying the 
ambiguities in summing the pcrturbative series one learns about the parametric dependence of 
power corrections on the hard scales and about their potential size[36]. 

By taking these points into account, it may be argued that the analysis of the higher- twist 
effects on the dependence of the meson wave function in single pseudoscalar and vector meson 
production at photon-photon collisions by the running coupling (RC) approach are significant 
from both theoretical and experimental points of view. 

In this work we will apply the running coupling approach[37] in order to compute effects of 
the infrared renormalons on the meson production in photon-photon collisions. This approach 
was employed also in our work[38] for calculation of the single meson production in proton- 
proton collisions. 

Photon-photon collisions represent a very useful tool for the study of hadron production. 
Basically, the more attractive feature is the simple, clean initial state, involving only QED 
interactions, which allows one to concentrates on the final hadronic state. This way, in fact, 
some of the more clean tests for pQCD models were proposed [14]. It is well known that exclusive 
77 — >■ hadron processes can be studied in the e+e" colliders, particularly 7*7* processes, play a 
spesial role in QCD [39] , since their analysis is under much better control than the calculation 
of hadronic processes, which requre the input of non-perturbative hadronic structure functions 
or wave functions. 

A precise measurement of the inclusive charged meson production cross section at ^/s = 
183GeV and ^/s = 2096*6^ is important for the photon-photon collisions program at the 
International Linear Collider (ILC). The results of our calculations are based on the photon- 
photon coUisions at ^/s — ISSGeV and a/s — 2096*6^. 

The higher-twist contributions to high-pT inclusive meson production in two-photon colli- 
sions, a single meson inclusive photoproduction and jet photoproduction cross sections were 
studied by various authors [40-42]. As experiments examining high-p^ particle production in 
two-photon collisions are improved, it becomes important to reassess the various contribu- 
tions which arise in quantum chromodynamics. Predicting for the higher-twist contributions, 
originally obtained in Ref.43, may now be refined using the exclusive-process QCD formalism 
devoloped in [44] . Another important aspect of this study is the choice of the meson model wave 
functions. In this respect, the contribution of the higher- twist Feynman diagrams to a single 
meson production cross section in photon-photon collisions has been computed by using various 
meson wave functions. Also, the leading and higher-twist contributions have been estimated 
and compared to each other. Within this context, this paper is organized as follows: in Sec. 



Hn we provide some formulas for the calculation of the contribution of the high twist diagrams. 
In Sec. mil we present formulas and an analysis of the higher-twist effects on the dependence of 
the meson wave function by the running coupling constant approach. In Sec. IIVI we provide 
the formulas for the calculation of the contribution of the leading-twist diagrams and in Sec. |Vl 
we present the numerical results for the cross section and discuss the dependence of the cross 
section on the meson wave functions. We state our conclusions in section |Vll 



II. CONTRIBUTION OF THE HIGH TWIST DIAGRAMS 

The higher-twist Feynman diagrams, which describe the subprocess 7g — )■ Mq contributes 
to 77 — )■ MX for the meson production in the photon-photon collision are shown in Fig. 1(a). 
The amplitude for this subprocess can be found by means of the Brodsky-Lepage formula [45] 

M(s,t)= / dxi I dx25il - xi - X2)^Mixi,X2,Q'^)TH{s,i;xi,X2). (2.1) 
Jo Jo 

In Eq.(2.1), Th is the sum of the graphs contributing to the hard-scattering part of the 
subprocess. The hard-scattering part for the subprocess under consideration is •yqi — > {qiq2)(l2, 
where a quark and antiquark form a pseudoscalar, color-singlet state {qiq2)- Here $(xi, X2, Q^) 
is the meson wave function, i.e., the probability amplitude for finding the valence qiq2 Fock 
state in the meson carry fractions Xi and X2-, Xi + X2 = 1. Remarkably, this factorization is 
gauge invariant and only requires that the momentum transfers in be large compared to 
the intrinsic mass scales of QCD. Since the distribution amplitude and the hard scattering 
amplitude are defined without reference to the perturbation theory, the factorization is valid to 
leading order in 1/Q, independent of the convergence of perturbative expansions. The Hard- 
scattering amplitude Th can be calculated in perturbation theory and represented as a series 
in the QCD running coupling constant asiQ"^)- 

The ^1^2 spin state used in computing may be written in the form 



V2 

PLheUcttyO, (2.2) 
^ T^, PT heHcity ± 1, 

where 6± = t(1/V2)(0, 1, ±i, 0) in a frame with {pm) 1,2 = and the N^^^^ project out a 
state of spins s, and pm is the four-momentum of the final meson. In our calculation, we have 
neglected the meson mass. Turning to extracting the contributions of the higher-twist subpro- 
cesses, there are many kinds of leading-twist subprocesses in 77 collisions as the background 



of the higher- twist subprocess •yq — )■ Mq, such as 7 + 7 — )■ q + q- The contributions from 
these leading-twist subprocesses strongly depend on some phcnomcnological factors, for exam- 
ple, quark and gluon distribution functions in meson and fragmentation functions of various 
constituents etc. Most of these factors have not been well determined, neither theoretically 
nor experimentally. Thus they cause very large uncertainty in the computation of the cross 
section of process 77 — > MX. In general, the magnitude of this uncertainty is much larger than 
the sum of all the higher-twist contributions, so it is very difficult to extract the higher-twist 
contributions. 

The Mandelstam invariant variables for subprocesses jq Mq are defined as 



(2.3) 



We have aimed to calculate the meson production cross section and to fix the differences 
due to the use of various meson model functions. We have used seven different wave functions: 
the asymptotic wave function (ASY), the Chernyak-Zhitnitsky wave function [6,9], the wave 
function in which two nontrivial Gegenbauer coefficients 02 and 04 have been extracted from 
the CLEO data on the 77* — > 7r° transition form factor [46], the Braun-Filyanov pion wave 
functions [11] and the Bakulev-Mikhailov-Stefanis pion wave function[47]. It should be noted 
that the wave functions of pions also are developed in Refs. [48-50] by the Dubna group. For p- 
meson wave function we used the Ball-Braun wave function[51]. 



Cl'\2x-l) + ^-Cr{2x-l) 



-,3/2. 



Cl'^{2x - 1) + 0.18(0.2) JC2'"(2a; - 1) 



.3/2, 



^BMs{x,l4) = "^asyix) [c^o^\2x - 1) + 0.188C2'/'(2a; - 1) - 0.13Cf\2x - 1)] , 
^CLEoix, III) = ^asy(x) [c^^\2x - 1) + 0.27C|/'(2x - 1) - 0.22cf '(2x - 1)] , 
^BF{x,fil) = ^asy{x) [Co^\2x - 1) + 0.44C|/'(2x - 1) + 0.25C|/'(2x - 1)] , 



C'J\2x - 1) = 1, CY\2x - 1) = ^{5{2x - If - 1), 



-,3/2. 



Cf\2x - 1) = ^(21(2x - 1)* - 14(2x - If + 1). 
8 

where /jr = 0.923GeV, /^=0.141 GeV, /J=0.16 GeV are the pion and p mesons decay 
constants. Here, we have denoted by x = xi, the longitudinal fractional momentum carried by 
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(2.4) 



the quark within the meson. Then, X2 = I — x and xi — X2 = 2x — 1. The pion and p meson 
wave functions is symmetric under the replacement Xi — X2 ^ X2 — Xi. 

Several important nonperturbative tools have been developed which allow specific predic- 
tions for the hadronic wave functions directly from theory and experiments. The QCD sum-rule 
technique and lattice gauge theory provide constraints on the moments of the hadronic distribu- 
tion amplitude. However, the correct meson wave function is still an open problem in QCD. It 
is known that the meson wave function can be expanded over the eigenfunctions of the one-loop 
Brodsky-Lepage equation, i.e., in terms of the Gegenbauer polynomials {Cn'^{'2x — 1)}, 



l+Y,ar.{Q')C'J'{2x-l) 



n=2.. 



(2.5) 



In the present work, we take into account the evolution of the meson wave function on the 
factorization scale. The evolution of the wave function on the factorization scale is governed 
by the functions a„((5^), 

In the case tt meson. 



«s(/^o) 



7n//30 



(2.6) 



72^ 50 74 ^ 364 ^ 
/3o 81' /3o 405' 
In Eq.(2.6), {7n} are anomalous dimensions defined by the expression. 



7n = Cp 



n+1 



(n+l)(n + 2) 



In the case p meson. 



(2.7) 



2^ll(±) 



Qs(Q 



2^^(7P^-7o)/(2/3o) 



Here, 



^11 — 

In 



1 - 



n+1 ^ 



> + l)(^ + 2) —J 



±^8 

In 3 



n+1 ^ 



i=2 



(2.8) 



The constants an{p^) = are input parameters that form the shape of the wave functions 
and which can be extracted from experimental data or obtained from the nonperturbative 



QCD computations at the normalization point /Xg. The QCD coupling constant asiQ"^) at the 
one-loop approximation is given by the expression 



Here, A is the fundamental QCD scale parameter, /3o is the QCD beta function one-loop 

coefficient, 

/3o = 11 - ^n/. 

The higher-twist subprocess 75 — >■ Mq contributes to 77 — >■ MX through the diagram of 
Fig. 1(a). We now incorporate the higher-twist (HT) subprocess 75 — > Mq into the full inclusive 
cross section. In this subprocess 75 Mq, photon and the meson may be viewed as an effective 
current striking the incoming quark line. With this in mind, we write the complete cross section 
in formal analogy with deep-inelastic scattering, 

E-^{-f^ MX) = - V / dx6{s + i + u)sGq/^{x, -i)^{^q Mq) + it ^ u), (2.10) 

dp TT _ Jo dt 

qq 

Here Gg/^ is the per color distribution function for a quark in a photon. The subprocess 
cross section for tt, pi and px production 

i(,,^M,) J ^^If*-^-*)!^^ + ''-^'"^ (2.11) 



where 



D(s, u) — eiu / dxi 



as{Ql)'^M{xi,Ql] 

Xi{l - Xi) 



+ 628 / dxi 



Xi(l - Xi] 



2\-\ 



(2.12) 



where Ql = s/2, Q\ = — m/2, represents the momentum squared carried by the hard gluon 
in Fig. 1(a), 61(62) is the charge of ^'1(^2) ^-^^d Cp — ^■ 

In the running coupling method this cross section was found in Ref.[52]. For pseudoscalar 
and longitudinally polarized meson 
da"^{ei,e2) _ 327r^CFaE 



dt 



9s2 

S2 



s 

26162 



^[Ili - 2h{hs + hu)l + lit] - ^AKli - 2K^{K^u + K2s)l+ 



t 



t 



t sut 

for the transversely polarized vector meson. 



da"'^{ei,e2) GAn^CpaE -i 



dt 



9s4 



— [eirt/2 - 625X2]' 



t 

(2.13) 
(2.14) 



Here, 



m = f f 

Jo Jo 



dxidx25{l - xi- X2)as{(l - Xi)s)^m(x, Q"^) 



X2 



and 



K2{u) 



r f 

Jo 



dxidx25{l - xi- X2)as{{l - Xi)s)^m{x, Q"^) 



X1X2 



1 pi 



dxidx2S{l — xi — X2)as{—xiu)^M{x, Q"^) 



JO 
1 rl 



dxidx2S{l - xi- X2)as{-Xiu))^M{x, Q^) 



IQ Jo X1X2 

The full cross section for tt and pi production is given by 



d'^p s -\-u 



^Gq/^{x,-i) 



3 



1 1 " 



+ 



s + t 



1 1 



In (2.15), the subprocess invariants are 



(2.15) 



s = xs, u = xu, t = t, 



(2.16) 



t = --^{xr - xf) = -mT\/se ^, 

with xr = {x% + x^Y^"^. Here xp = 2(pm)||/-s/s and xt = '2{pm)±/Vs = Iprj \fs specify the 
longitudinal and transverse momentum of the meson. In terms of these the rapidity of M is 
given by 

y = + ^f)/{xr - Xf)] 

where rriT is the transverse mass of meson, which is given by 

As seen from (2.11) the subprocess cross section for longitudinal pL production is very similar 
to that for tt production, but the transverse pt subprocess cross section has a quite different 
form. 

Let as first consider the frozen coupling approach. In this approach we take equal the 
four-momentum square Q\ 2 of the hard gluon to the meson's transverse momentum square 
Q\,2 — Pt- this case the QCD coupling constant ag in the integral (2.12) does not depend on 



integration variable. After this substitution, calculation of integral (2.12) becomes easy. Hence, 
the effective cross section obtained after substitution of the integral (2.12) into the expression 
(2.15) is referred as the frozen coupling effective cross section. We will denote the higher-twist 
cross section obtained using the frozen coupling constant approximation by (E^^)'^. 

For a full discussion, we consider a difference A^^ between the higher- twist cross section 
combinations and 

A^^^ = S^^ - S^^_ = ^M.^(77 ^ M^X) - EM-^h7 ^ M'X). (2.17) 

We have extracted the following higher-twist subprocesses contributing to the two covariant 
cross sections in Eq.(2.11) 

IQi ^ {QiQ2)Q2 ,7Q2^ {QiQ2)Q2 (2-18) 

As seen from Eq.(2.15), at fixed pr, the cross section falls very slowly with s. Also, at fixed 
s, the cross section decreases as l/p^, multiplied by a slowly varying logarithmic function which 
vanishes at the phase-spase boundary. Thus, the pr spectrum is fairly independent of s expect 
near the kinematic limit. 



III. THE RUNNING COUPLING APPROACH AND HIGHER- TWIST MECHA- 
NISM 

In this section we shall calculate the integral (2.12) using the running coupling constant 
approach and also discuss the problem of normalization of the higher-twist process cross section 
in the context of the same approach. 

As is seen from (2.12), in general, one has to take into account not only the dependence of 
Oi{Qi^2) on the scale (5i,2) but also an evolution of ^{x, Qi,2) with 2- The meson wave function 
evolves in accordance with a Bethe-Salpeter-type equation. Therefore, it is worth noting that, 
the renormalization scale (argument of ag) should be equal to Qf — X2S, Q\ — —Xiil, whereas 
the factorization scale [Q^ in $M(a^, Q^)] is taken independent from x, we take = p^. 
Such approximation docs not considerably change the numerical results, but the phenomenon 
considered in this article (effect of infrared renormalons) becomes transparent. The main 
problem in our investigation is the calculation of the integral in (2.12) by the running coupling 
constant approach. The integral in Eq.(2.12) in the framework of the running coupling approach 
takes the form 



Jo 



1-x 



(3.1) 



The as{XiJ^) has the infrared singularity atx— >1, ifA = l — xorx— >0, ifA = a; and 
as a result integral (3.1) diverges (the pole associated with the denominator of the integrand 
is fictitious, because $m ~ (1 ~ 2;), and therefore, the singularity of the integrand at x = 1 is 
caused only by as((l — x)/x|,y)). For the regularization of the integral we express the running 
coupling at scaling variable asi^fij^^) with the aid of the renormalization group equation in 
terms of the fixed one as{Q'^). The renormalization group equation for the running coupling 
a = cKs/tt has the form [31] 



where 



dlnX - 4 ^"^^"^ 



3o = 11 - 



(3.2) 



The solution of Eq.(3.2), with the initial condition 



a(A)|A=i = a = as{Q^)/n, 



is [31] 



a{X) 



a 



1 + a—lnX 
4 



(3.3) 



This transcendental equation can be solved iteratively by keeping the leading aHn^X order. 
This term is given by 



oisiXQ^) 



asm 

1 + InX/t 



After substituting Eq.(3.4) into Eq.(2.12) we get 

'^^^a,( A/x|J$M(a:,(3') 
x(l — x) 



D(s, u) — tit I dx- 
Jo 

eitas\s) / ax—. ; — — ; — - -\- e2uas\—u) / ax—. —r- 

' x{l-x){l + lnX/h) Vo x{l-x){l 



x[\ — x) 

x,Q') 



+ lnX/t2) 



eita. 



e2uas 



.1 ^asy{x) 

,{s) I dx 

^0 

{—u) / dx 

^0 



1 + 5^2,4,.. (^nilJ'O 



cT{^x - 1) 



a;(l -a;)(l + /nA/ti) 



1 + Em,..^^(^o) 



"s(Mo) 



Cl'\2x - 1) 



x{l- x){l + lnX/t2) 



(3.4) 



(3.5) 



where ti = 47r/as(s)/3o, ^2 = ^7r/as{—u)l3o. The integral (3.5) is common and, of course, still 
divergent, but now it is recast into a form, which is suitable for calculation. Using the running 
coupling constant approach, this integral may be found as a perturbative series in as 

oo 

(3-6) 

n=l 

The expression coefficients Sn can be written as power series in the number of light quark 
flavors or, equivalently, as a series in power of Po, as Sn — Cn/^o"^- The coefficients C„ of this 
series demonstrate factorial growth Cn ~ {n — l)\, which might indicate an infrared renormalon 
nature of divergences in the integral (3.5) and corresponding series (3.6). The procedure for 
dealing with such ill-defined series is well known; one has to perform the Borel transform of the 
series [53] 

oo 



nl 

n=0 



then invert B[D]{u) to obtain the resummed expression (the Borel sum) D{s,u). After this 
we can find directly the resummed expression for D{s,u). The change of the variable x to 
z — ln{l — x), as ln{l — x) — InX. Then, 

D{s,u) ^ eitas{s)ti — --— — - + e2uas{-u)t2 — ry— — r (3.7) 

Jo X{1 - X){ti + Z) Jo X{1 - X){t2 + Z) 

For the calculation the expression (3.7) we will apply the integral representation of l/(t-|-^) 
[54]. After this operation, formula (3.7) is simplified and we can extract the Borel sum of the 
perturbative series (3.6) and the corresponding Borel transform in dependence from the wave 
functions of the meson, respectively. Also after such manipulations the obtained expression can 
be used for numerical computations. 

It is convenient to use the following integral representation for l/{t + z): 

/ e-^'+'^du (3.8) 
t + z Jo ^ ^ 

After inserting Eq.(3.8) into (3.7), then, we obtain 

n.-^ f r^,, r ^Mix,Q')e'^'^+^>dudx , 

D{s,u) = eitas{s)ti / / r h 

Jo Jo x[l — X) 

. . f' r ^M{x,Q')e-('^+^^^dudx 
e2uas[-u)t2 / / r . (3.9) 



Jo 



x{l — x) 



In the case of ^asy{x) for /i(s), /2(s), Ki{u), K2{u), we get 

1 



1 — u 2 — u 



his) 



Ki{u) 



Jo 



due 



-tiu 



1-U 



12 



Jo 



due 



-t2U 



due 



1 — u 2 — u 
1 



-t2U 



/3o Jo 

In the case of the Q^) wave function, we find 



1-u 



0.27(0.3) 



(50/162), (26/162) 



due-^'"" 
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1 — u 2 — u 



24 40 

+ 



20 



0.27(0.3) 



m - 



«s()"o) 



(50/162), (26/162) 



1 — M 2 — M 3 — ti 4 — ti_ 
1 



'0 
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1-u 
20 
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20 



0.27(0.3) 



K,{U) : 

«.(Q^) ' 



4v^7r/, 



•0 Jo 



due 



1 — u 2 — u 3 — u 
1 1 



-t2U 



(50/162), (26/162) 



1 — u 2 — u 

24 40 

+ 



20 



0.27(0.3) 



K2{u) 
."s(/^o) 



4y67r/p 
Po Jo 

(50/162), (26/162) 



1 — u 2 — u 3 — u A — u_ 
1 



due-^^"" 
4 



1 -M 

20 



20 



1 — u 2 — u 3 — M 



In the case of the ^cleo{x, Q^) wave function, we get 



Po 7o 



1 



24 40 

+ 
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Equation(3.1) and (3.2) is nothing more than the Borel sum of the perturbative series (3.6), 
and the corresponding Borel transform in the case ^asy{x) is 
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The series (3.6) can be recovered by means of the following formula: 

n-l 



The Borel transform B[D]{u) has poles on the real u axis &t u — 1; 2; 3; 4; 5; 6, which confirms 
our conclusion concerning the infrared renormalon nature of divergences in (3.6). To remove 
them from Eqs. (3. 10-3.20) some regularization methods have to be applied. In this article we 
adopt the principal value prescription. We obtain: in the case ^asy 
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where Li (A) is the logarithmic integral for A > 1 defined as the principal value [55] 
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Hence, the effective cross section obtained after substitution of the expressions (3.10-3.13) 
into the expression (2.15) is referred as the running coupling effective cross section. We will 
denote the higher-twist cross section obtained using the running coupling constant approach 
by (Ef^^)'--. 



IV. CONTRIBUTION OF THE LEADING-TWIST DIAGRAMS 



Regarding the higher-twist corrections to the meson production cross section, a comparison 
of our results with leading-twist contributions is crucial. The contribution from the leading- 
twist subprocess 77 — )■ qq is shown in Fig. 1(b). The corresponding inclusive cross section for 
production of a meson M is given by 



da 



-E f^S{s + i + u)sDfiz,-iAjj^qq) (4.1) 

jj^MX ^ g^q Jo Z dt 



where 



t ^ u 
s — s, t — - u — — 

z z 



Here s, t, and u refer to the overall 77 MX reaction. D^[z, —t) represents the quark 
fragmentation function into a meson containing a quark of the same flavor. For 7r+ production 
we assume /„ = In the leading- twist subprocess, meson is indirectly emitted from 

the quark with fractional momentum z. The 5 function may be expressed in terms of the parton 
kinematic variables, and the z integration may then be done. The final form for the leading- 
twist contribution to the large-pr meson production cross section in the process 77 — > MX 
is 
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We should note that D{z, —t)/z behaves as 1/2;^ as 2; ^ 0. For the kinematic range considered 
in our numerical calculations, D[z,—t)/z increases even more rapidly. We obtain from the 
final cross section Eq.(4.2), following conclusion: At fixed pt, the cross section decreases with 
s asymptotically as 1/s. At fixed s, the D{z, —i) function causes the cross section to decrease 
rapidly as pr increases towards the phase-spase boundary (2; — > 1). As s increases, the phase- 
spase boundary moves to higher px, and the px distribution broadens. 



V. NUMERICAL RESULTS AND DISCUSSION 



In this section, we discuss the numerical results for higher-twist effects with higher-twist 
contributions calculated in the context of the running coupling constant and frozen coupling 



approaches on the dependence of the chosen meson wave functions in the process 77 MX. 
We have calculated the dependence on the meson wave functions for the higher-twist contribu- 
tion to the large-pT single pseudoscalar tt"*" and vector p^, mesons production cross section 
in the photon-photon collision. The 7r~, cross sections are, of course, identical. In 

the calculations, we use the asymptotic wave function ^asy, the Chernyak-Zhitnitsky ^cz, the 
pion wave function from which two nontrivial Gegenbauer coefficients 02 and 04 have been 
extracted from the CLEO data on the 7r°7 transition form factor [46], the Braun-Filyanov pion 
wave functions [11], and the Bakulev-Mikhailov-Stefanis pion wave function[ 47]. For p-meson 
we used Ball-Braun wave function[51]. For the higher-twist subprocess, we take •jqi — )■ (^'1^2)52, 
IQ2 ~^ {0.10.2)0.2 contributing to 77 — >■ MX cross sections. Inclusive meson photoproduction 
represents a significant test case in which higher-twist terms dominate those of leading-twist 
in certain kinematic domains. For the dominant leading-twist subprocess for the meson pro- 
duction, we take the photon-photon annihilation 77 qq, in which the M meson is indirectly 
emitted from the quark. For example, the quark distribution function inside the photon has 
been used [56]. The higher-twist subprocesses probe the meson wave functions over a large 
range of squared momentum transfer, carried by the gluon. Therefore, in the diagram given 
in Fig. la we take Ql — X2S, Q\ — —Xiu , which we have obtained directly from the higher- 
twist subprocesses diagrams. The same 2 has been used as an argument of ols{Q\^-^ in the 
calculation of diagram. 

The results of our numerical calculations are plotted in Figs. 2-33. First of all, it is very 
interesting to compare the resummed higher-twist cross sections with the ones obtained in the 
framework of the frozen coupling approach. In Figs. 2-4 we show the dependence of higher- 
twist cross sections (E^+)° calculated in the context of the frozen couphng, (E^+)''^* in the 
context of the running coupling constant approaches and also the ratio it! = (E^+)''^*/E^+)° 
as a function of the meson transverse momentum pt for different meson wave functions at 
1/ = 0. It is seen that the values of cross sections (E^+)°, (S^+)*'^*, and R for fixed y and 
\/s depend on the choice of the meson wave function. As seen from Figs. 2-3 in both cases, 
frozen coupling and running coupling constant approaches the higher-twist differential cross 
section is monotically decreasing with an increase in the transverse momentum of the meson. 
As is seen from Fig.4, when the transverse momentum of the meson is increasing, the ratio R is 
decreasing. But, as shown in Fig.4, in the region 5 GeV/c < < 80 GeV/c higher-twist cross 
section calculated in the context of the running coupling method is suppressed by about 2-4 
orders of magnitude relative to the higher-twist cross section calculated in the framework of the 
frozen coupling method. In Figs. 5 and 6, we shows the dependence of the ratio (E^+)°/(E^) 



and (S^+)^'^*/(S^^) as a function of the meson transverse momentum px for different meson 
wave functions. Here (S^^) is the leading- twist cross section, respectively. As seen from 
Fig.6, in the region 5 GeV/c < pt < 10 GeV/c higher-twist cross section for $_BM5(a^, <5^)) 
and ^cLEoix^Q"^)) wave functions calculated in the context of the running coupling method 
is suppressed by about one orders of magnitude relative to the leading-twist cross section, 
but in the region 10 GeV/c < < 90 GeV/c ratio is decreasing with an increase in the 
transverse momentum of the meson. In Figs. 7-9 we show the dependence (A^)^, {A^^'^Y^'^, 
and the ratio r=(A^^)''^*/(A^"^)°, as a function of the meson transverse momentum pt for the 
different meson wave functions. Here, (Af/)0 = {^f/+Y - (Sm-)° and (Af/)'''^^ = {^fjlY"' - 
(S^^)^^'*. As seen from Figs. 7 and 8, the difference of the higher-twist differential cross section 
is decreasing with an increase in the transverse momentum of the meson. The dependence, as 
shown in Fig.9, is identically equivalent to Fig.4. In Figs. 10-17, we have depicted higher-twist 
cross sections, ratios (Sf,^)", R = (S^,^ )-7(E^^ )«, r=(A^^)-V(Af/)0, (Af/)", 

(Af/)'^'=^ (S^^^)7(S^j;) and (Sf^'^)'-^V(S^7+) as a function of the rapidity y of the meson at 
^/s = 183 GeV and pr = 14.6 GeV/c. At = 183 GeV and pr = 14.6 GeV/c, the meson 
rapidity lies in the region —2.52 < y < 2.52. 

As seen from Fig. 10 and Fig.l4, in the region (—2.52 < y < —1.92), the cross section for all 
wave functions increases with an increase of the y rapidity of the meson and have a maximum 
approximately at the point y = —1.92. Besides that, the cross sections decrease with an increase 
in the y rapidity of the meson. But, as seen from Figs. 12 and 13 in the region (—2.52 < y < 1.92) 
the cross section for all wave functions increases with an increase of the y rapidity of the meson 
and has a maximum approximately at the point y = 1.92. But, as seen from Figs. 16-17 in the 
region (—2.52 < y < 1.92) the cross section for all wave functions has a minimum approximately 
at the point y = 1.92. As is seen from Figs. 10-17, cross sections, the ratios R and r are very 
sensitive to the choice of the meson wave functions. It should be noted that the magnitude 
of the higher-twist cross section for the pion wave functions ^BMsi^jQ"^) and ^cLEoi^jQ"^) is 
very close to the asymptotic wave function ^asy{x)- Also, the distinction between R{^asy{x)) 
with R{<!>BMs{x, Q')), R{<l>cz{x, Q2)), R{<l>cLEo{x, Q')), R{<I>bf{x, Q2)), R{<l>BBiL){x, Q^)) and 
R{^BB{T){x, Q"^)) have been calculated. For example, in the case of a/s = 183 GeV , y = 0, the 
distinction between R{<!>asy{x)) with R{<l>i{x,Q^)) [i=BMS, CZ, CLEO, BF, BB(L), BB(T)] as 
a function of the meson transverse momentum pt is shown in Table [B Thus, the distinction 
between R{^asyix)) with R{^i{x,Q'^)), [i=BMS, CLEO] is maximum at pr = 5 GeV/c, with 
R{^cz{x)) at Pt = 50 GeV/ c; the distinction between R{^asy{x)) with R{^bf{x,Q^)), is 
maximum at pt = 90 GeV/c; but the distinction R{^asy{x)) with i?($j(x, Q^)), [i=BB(L), 



BB(T)] is maximum at px = 75 GeV/c, Also, we have calculated the distinction between 
R{^asy{x)) with R{^i{x,Q^)) [i=BMS, CZ, CLEO, BF, BB(L), BB(T)] as a function of the 
rapidity y of the meson. For example, in the case of ^/s = 183GeV, Pt = 14:.QGeV/c the 
distinction between R{^asy{x)) with R{^i{x,Q^)) [i=BMS, CZ, CLEO, BF, BB(L), BB(T)] 
as a function of the rapidity y of the meson is presented in Table [III 

We have also carried out comparative calculations in the center-of-mass energy ^/s = 
209 GeV. The results of our numerical calculations in the center-of-mass energies a/s = 
209 GeV are plotted in Figs. 18-33. Analysis of our calculations at the center-of-mass ener- 
gies a/s = 183 GeV and y/s = 209 GeV, show that with the increase in beam energy val- 
ues of the cross sections, ratio R = (S^+)^^*/(S;^+)°, and contributions of higher-twist to 
the cross section decrease by about 1-2 order. Therefore the experimental investigation of 
higher-twist effects include renormalon effects conveniently in low energy. On the other hand, 
the higher-twist corrections and ratios R and r are very sensitive to the choice of the meson 
wave function. Also, the distinction between R{^asy{x)) with R{^bms{x, Q^)), R{^cz{x, Q^)), 
Ri^CLEoix,Q^)), R{^BFix,Q^)), Ri(^BB{L)ix,Q^)) and Ri(^BB{T)ix,Q^)) have been calcu- 
lated. For example, in the case of a/s = 209 GeV, y = 0, the distinction between R{^asy{x)) 
with R{^i{x,Q^)) [i=BMS, CZ, CLEO, BF, BB(L), BB(T)] as a function of the meson 
transverse momentum px is shown in Table IIIII Thus, the distinction between R{^asy{x)) 
with R{^i{x,Q'^)), (i=BMS, CLEO) is maximum at px = 10 GeV/ c, with R{^czix)) at 
Pt = 65 GeV/ c; the distinction between R{^asy{x)) with R{^bf{x,Q'^)), is maximum at 
Pt = 100 GeV/c; but the distinction R{^asy{x)) with R{^i{x,Q^)), (i=BB(L), BB(T)) is 
maximum at px = 65 GeV/c. Also, we have calculated the distinction between R{^asy{x)) 
with i?($i(x,g2)) [i=BMS, CZ, CLEO, BF, BB(L), BB(T)] as a function of the rapidity y 
of the meson. For example, in the case of ^/s = 209GeV, pt = lQ.7GeV/c the distinction 
between R{^asy{x)) with R{^i{x,Q^)) [i=BMS, CZ, CLEO, BF, BB(L), BB(T)] as a func- 
tion of the rapidity y of the meson is presented in Table IIVI The calculations show that the 
ratio R{^i{x,Q^))/R{^asy{x)), (i=CLEO, CZ, BF, BMS, BB(L), BB(T)) for all values of the 
transverse momentum pt of the meson identically equivalent to ratio r{^i{x,Q'^))/r{^asy{x))- 
Results of our numerical calculations demonstrate that in the renormalon approach there are 
not difference between results obtained using the cross sections (2.13), (2.14) or (2.11), (2.12). 

The total integrated luminosity of LEP is 612.8p6~^ and total luminosity of ILC required 
is 500/6~^, also a peak luminosity of ILC is 1000 fb~^ during the first phase of operation at 
209GeV -^500GeV . In our calculations of the higher-twist cross section of the process the depen- 
dence of the transverse momentum of meson appears in the range of (10~^° -^ 10~^^)mb/ GeV"^ , 



or (10 ^ -^ 10 ^^)pb/GeV'^. Therefore, higher-twist cross section obtained in our paper should 
be observable at LEP and ILC. 



VI. CONCLUDING REMARKS 



In this work we have calculated the single meson inclusive production via higher-twist mech- 
anism and obtained the expressions for the subprocess — > Mq cross section for mesons with 
symmetric wave functions. For calculation of the cross section we have applied the running cou- 
pling constant method and revealed infrared renormalon poles in the cross section expression. 
Infrared renormalon induced divergences have been regularized by means of the principal value 
prescription and the resummed expression (the Borel sum) for the higher-twist cross section has 
been found. The higher-twist cross sections were calculated in the frozen coupling and running 
coupling approaches. The resummed higher-twist cross section differs from that found using the 
frozen coupling approach, in some regions, considerably. Also we demonstrated that higher- 
twist contributions to single meson production cross section in the photon-photon collisions 
have important phenomenological consequences. We have obtained very interesting results. 
The ratio R for all values of the transverse momentum pr and of the rapidity y of the meson 
identically equivalent to ratio r. Our investigation enables us to conclude that the higher- twist 
meson production cross section in the photon-photon collisions depends on the form of the 
meson model wave functions and may be used for their study. Analysis of our calculations 
shows that the magnitude of cross sections of the leading-twist is larger than the higher-twist 
cross sections ones calculated in the frozen coupling approach in 2-4 order. But, in some re- 
gions of transverse momentum of the meson, the higher-twist cross section calculated in the 
context of the running coupling method is comparable with the cross sections of leading- twist. 
Further investigations are needed in order to clarify the role of high twist effects in this process. 
We have demonstrated that the resummed result depends on the meson model wave functions 
used in calculations. The production of high-p^ meson probes the short-distance dynamics of 
photon-photon reactions. In addition to providing tests of perturbative QCD, 77 processes 
with real or almost real photons give us information on the photon structure function which is 
complementary to the information gained from deep inelastic scattering on a real photon. The 
latter process essentially probes the quark distribution while high-p^ meson production is also 
sensitive to the gluon distribution of the photon. As it is well known high-p^ processes induced 
by real photons have a rather complex structure. This arises from the fact that the photon 
couples to the hard subprocess either directly or through its quark and gluon content. In par- 



ticular, meson production in photon-photon coUisions takes into account infrared renormalon 

effects: this opens a window toward new types of photon structure function which can not be 
measured by the lepton-photon scatterings. 
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TABLE I: The distinction between R{^asy{x)) with R{^i{x,Q'^)) (i=BMS, CZ, CLEO, BF, BB(L), 
BB(T) ) at cm. energy = 183GeF. 
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TABLE IL The distinction between R{<i>asy{x)) with R{^i{x,Q'^)) (i=BMS, CZ, CLEO, BF, BB(L), 
BB(T)) at cm. energy = 183Gey. 
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TABLE III: The distinction between R{^asy{x)) with R{^i{x,Q'^)) (i=BMS, CZ, CLEO, BF, BB(L), 
BB(T)) at cm. energy = 209GeF. 
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TABLE IV: The distinction between R{^asy{x)) with R{^i{x,Q'^)) (i=BMS, CZ, CLEO, BF, BB(L), 
BB(T)) at cm. energy = 209GeF. 
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FIG. 1: (a): The higher-twist contribution to 77 — >• MX; (b): The leading-twist contribution to 

77 MX 
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FIG. 2: Higher-twist M production cross section (S^-^)'^ as a function of thepT transverse momentum 
of the meson at the cm. energy ^/s = 183 GeV . 
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FIG. 3: Higher-twist M production cross section {Y^HTyes ^ ^ function of the pT transverse momen- 
tum of the meson at the cm. energy ^/s = 183 GeV . 
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FIG. 4: Ratio R = (Sf/)''^7(Sf/)°, where higher-twist contribution are calculated for the meson 
rapidity y = at the cm. energy ^/s = 183 GeV as a function of the meson transverse momentum, 
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FIG. 5: Ratio (S^ as a function of the pT transverse momentum of the meson at the cm. 

energy -y^s = 183 GeV. 
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FIG. 6: Ratio {Ti^fY^'^/ as a function of the pT transverse momentum of the meson at the cm. 
energy =183 GeV . 
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FIG. 7: The difference of the higher-twist cross section, (A|^^)° = - as a function 

of the meson transverse momentum, px, at the cm. energy y/s = 183 GeV . 
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FIG. 8: The difference of tlie liigher-twist cross section, (A^"^'*''^* 



/'V-n-i \res /v 

M J — \^M+) ~ y^M 
function of the meson transverse momentum, px, at the cm. energy ^/s = 183 GeV. 
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FIG. 9: Ratio r = (A|P')'''^*/(A^^)'', where higher-twist contributions are calculated for the meson 
rapidity y = at the cm. energy \/s = 183 GeV , as a function of the meson transverse momentum, 
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FIG. 10: Higher-twist M production cross section (S^-^)", as a function of the y rapidity of the meson 
at the transverse momentum of the meson px = 14.6 GeV/c, at the cm. energy y/s = 183 GeV. 



> 
O 



10"' 



10" 



10"^*^ b- 



X s 



10" 



10"" 




BMS 

cz 

CLEO p^=14.6GeV/c 
BF 
BB(L) 
BB(T) 



-1 





y 



FIG. 11: Higher-twist M production cross section as a function of the y rapidity of the meson 

at the transverse momentum of the meson px = 14.6 GeV/c, at the cm. energy y/s = 183 GeV. 



10° 



10^ r 



X a 



10' - 



X a 



10" r 



10"^ 



-3 






y 



FIG. 12: Ratio R = {Tj^^Y'^^ / {Y^^^)^ , as a function of the y rapidity of the meson at the transverse 
momentum of the meson pT = 14.6 GeV/c, at the cm. energy y/s = 183 GeV. 



10" w 



10^ r 



3= a 
< 



10^ t 

£ a 
<1 



10" r 



10"' 




/ y 



-2 





y 



FIG. 13: Ratio r = (A^-^)'"'^*/(A^-^)'^, as a function of the y rapidity of the meson at the transverse 
momentum of the meson = 14.6 GeV/c, at the cm. energy = 183 GeV 




FIG. 14: The difference of the higher-twist cross section, (AI^'^)" = — (E^T!)*^, as a function 

of the y rapidity of the meson at the transverse momentum of the meson px = 14.6 GeV/c, at the 
cm. energy ^/s = 183 GeV . 
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FIG. 15: The difference of the higher-twist cross section, (A 



HT\res 
M ) 



HT \res 
M+) 



HT \res 
M-> ' 



as a 



function of the y rapidity of the meson at the transverse momentum of the meson px = 14.6 GeV jc^ 
at the cm. energy = 183 GeV . 




FIG. 16: Ratio {Jl^f)^/ (S^J), as a function of the y rapidity of the meson at the transverse momentum 
of the meson pT = 14.6 GeV/c, at the cm. energy y/s = 183 GeV . 
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FIG. 17: Ratio {Jlj^jY^^/iTi^), as a function of the y rapidity of the meson at the transverse 
momentum of the meson px = 14.6 GeV/c, at the cm. energy y/s = 183 GeV . 



-IQ-^" I 1 I 1 I 1 I 1 I 1 I 

20 40 60 80 100 

, GeV/c 

FIG. 18: Higher-twist M production cross section as a function of the pt transverse momentum 

of the meson at the cm. energy y/s = 209 GeV . 
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FIG. 19: Higher-twist M production cross section {Jlj^jY^^ as a function of the pT transverse mo- 
mentum of the meson at the cm. energy ^/s = 209 GeV . 
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FIG. 20: Ratio R = (Sf/)'^'=7(Sf/)°, where higher-twist contribution are calculated for the meson 
rapidity y = at the cm. energy ^/s = 209 GeV as a function of the meson transverse momentum, 
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FIG. 21: Ratio )'^/ (S^j ), as a function of the pT transverse momentum of the meson at the cm. 

energy y/s = 209 GeV . 
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FIG. 22: Ratio as a function of the pT transverse momentum of the meson at the 

cm. energy ^/s = 209 GeV .. 
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FIG. 23: The difference of the higher-twist cross section, (A|^-^)'' = — as a function 

of the meson transverse momentum, pT-, at the cm. energy a/s = 209 GeV . 
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FIG. 24: The difference of the higher-twist cross section, {A^^Y^"^ = {T.^l+Y"' - {Y.^_Y''% 
as a function of the meson transverse momentum, px, at the cm. energy = 209 GeV. 
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FIG. 25: Ratio r = (Af/ )''^7(Af/)°, where higher-twist contributions are calculated for the meson 
rapidity y = at the cm. energy ^/s = 209 GeV, as a function of the meson transverse momentum, 
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FIG. 26: Higher-twist M production cross section (S^-^)", as a function of the y rapidity of the meson 
at the transverse momentum of the meson pT = 16.7 GeV/c, at the cm. energy ^/s = 209 GeV. 
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FIG. 27: Higher-twist M production cross section {T^^.j'^Y^'^ , as a function of the y rapidity of the meson 
at the transverse momentum of the meson px = 16.7 GeV/c, at the cm. energy ^/s = 209 GeV. 




FIG. 28: The difference of the higher-twist cross section, (A^^-^)'' = — (S^-^)*^, as a function 

of the y rapidity of the meson at the transverse momentum of the meson px = 16.7 GeV/c, at the 
cm. energy ^/s = 183 GeV . 



10" 



> 

a 



10"' 



E 1 


' 1 


' 1 ' 1 ' : 


















' if/ 












: ■'/I./ / 

: //■/ ' 

■ 'I / 


asy 




BMS 

CZ 




; 


CLEO 


P^=16.7GeV/c ; 


_ ; 


BF 




; ( 


BB(L) 




_ / 


BB(T) 




1 


1 






-2 



as a 



FIG. 29: The difference of the higher-twist cross section, (A^^-^)^'^'' = 
function of the y rapidity of the meson at the transverse momentum of the meson px = 16.7 GeV/c. 
at the cm. energy ^/s = 209 GeV . 
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FIG. 30: Ratio R = {Ti^^Y'^'^ / {Ti^^)^ , as a function of the y rapidity of the meson at the transverse 
momentum of the meson pT = 16.7 GeV/c, at the cm. energy ^/s = 209 GeV. 
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FIG. 31: Ratio {Ti^f)^/ (S^J^), as a function of the y rapidity of the meson at the transverse momentum 
of the meson px = 16.7 GeV jc, at the cm. energy ^/s = 209 GeV. 
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FIG. 32: Ratio {T^^f)^/ (S^J"), as a function of the y rapidity of the meson at the transverse momentum 
of the meson px = 16.7 GeV/c, at the cm. energy y/s = 209 GeV 




FIG. 33: Ratio r = {Af.jY^^/{A^,j)^, as a function of the y rapidity of the meson at the transverse 
momentum of the meson pT = 16.7 GeV/c, at the cm. energy y/s = 209 GeV 



